Model-dependent high-energy neutrino flux from Gamma-Ray Bursts 
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The IceCube Collaboration recently reported a stringent upper limit on the high energy neutrino 
flux from GRBs, which provides a meaningful constraint on the standard internal shock model. 
Recent broad band electromagnetic observations of GRBs also challenge the internal shock paradigm 
for GRBs, and some competing models for 7-ray prompt emission have been proposed. We describe 
a general scheme for calculating the GRB neutrino flux, and compare the predicted neutrino flux 
levels for different models. We point out that the current neutrino flux upper limit already disfavors 
any dissipative photosphere model that invokes high energy proton acceleration at the same site, 
and challenges the internal shock model. If the neutrino flux upper limit continues to go down in 
the next few years, then it would suggest that the GRB emission site is at a larger radius than the 
internal shock radius or for some reason protons do not get accelerated to high energies at the site 
where 7-ray photons are produced. The larger-radius solution for the low neutrino flux from GRBs 
would provide support to magnetic dissipation models that invoke a large dissipation radius, such 
as the ICMART model. 

PACS numbers: 95.55.Vj; 95.85.Ry; 98.70.Rz 



I. Introduction. — As energetic, non-thermal photon 
emitters, gamma-ray bursts (GRBs) have long been re- 
garded as efficient cosmic ray accelerators [l| . If protons 
in a GRB jet can be accelerated to an energy E p so that 
the condition 
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(1) 

is satisfied, significant neutrino emission is possible via 
the j>7 mechanism at the A-resonance. Here T is the bulk 
Lorentz factor, E 1 is photon energy in observer frame, 
toa = 1-232 GeV and m p = 0.938 GeV are the rest 
masses of A + and proton, respectively. For GRBs, a 
guaranteed target photon source for j>7 interaction is the 
burst itself. For the typical peak photon energy E 1 ~ 
several hundred keV, the corresponding neutrino energy 



E v ~ 0.05£„ 



(2) 



is in the sub-PeV regime Q which is well suited for detec- 
tion with the current generation of high energy neutrino 
detectors such as the IceCube Q. Indeed over the years, 
the IceCube Collaboration have been searching for high 
energy neutrino signals coincident with GRBs in time and 
direction, and progressively deeper non-detection upper 
limits have been placed [H, H[ , which are now beginning 
to constrain the standard GRB internal shock model Q ■ 
The current IceCube upper limit was claimed to be at 
least a factor of 3.7 smaller than the theoretical pre- 
dictions for neutrino flux from GRBs according to the 
internal-shock model, which has raised further doubt re- 
garding the viability of GRBs as sources of UHECRs @ . 
More detailed, follow-up, calculations 041] suggest that 
the current limit is still not deep enough to provide sig- 
nificant constraints on the validity of the internal shock 
model. However, the model would be severely challenged 
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if the upper limit continues to go down in the next few 
years. 

On the other hand, besides the internal shock model, 
there have been other GRB prompt emission models pro- 
posed oyer the years, including dissipative photosphere 
11| and magnetic dissipation models at large 
13j. Recent GRB observations with Swift and 
Fermi missions have suggested that the simple fireball 
internal shock model cannot account for the rich data 
of GRB prompt emission [14| . and these other mecha- 
nisms for GRB prompt emission become more attractive. 
The neutrino signal predictions of these prompt emission 
models could be very different from what is predicted 
for the internal shock model. The progressively strin- 
gent upper limit of neutrino flux would start to rule out 
some of these models. Since different models invoke dif- 
ferent assumptions regarding the GRB jet composition, 
future neutrino data would hold the key to diagnose the 
poorly known GRB jet composition. In this paper, we de- 
velop a general method for calculating the neutrino flux 
for a wide variety of GRB prompt emission models, and 
discuss how the current upper limit constrains various 
models. 

77. General formalism. Our general formalism closely 
follows the notations adopted by the IceCube Collabo- 
ration [4j, but we make the following changes: (1) In all 
the previous GRB neutrino flux calculations, the internal 
shock model has been implicitly assumed, so that the ra- 
dius where protons are accelerated and the radius where 
7-ray photons are generated are both taken as 



R = Ris =T 2 c6t min /(l + z), 



(3) 



where St mm is the minimum variability time scale ob- 
served in a GRB light curve. This expression, even 
though widely used, is not relevant for models other 
than the internal-shock model. For instance, in the 
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dissipative photosphere models, the photosphere radius 
R p h < i?is, and St mm could reflect the intrinsic vari- 
ability time scale of the central engine, which could 
be larger than the angular spreading time defined by 
-Rph/(r 2 c). The magnetic dissipation models (e.g. the 
Internal Collision-induced MAgnetic Reconnection and 
Turbulence or ICMART model P3) can have a GRB 
emission site R > R\s- The rapid variability time scale 
6t in these models is related to the time scale of relativis- 
tic mini-jets in the emission reg ion driven by relativistic 
turbulence or reconnection [12| . [iH . [TEj . To account for 
these possibilities, in our formalism we consider the pri- 
mary parameters to be R and T instead of St and T. (2) 
In the internal shock model, 7-rays and neutrinos are gen- 
erated by electrons and protons accelerated by the same 
shocks. And a parameter f e (electron-to-proton energy 
ratio in the internal shocks) relates the neutrino flux to 
the observed 7-ray flux - the assumption being that elec- 
tron energy is efficiently converted to 7-rays and neutrino 
flux is proportional to the energy in high energy protons. 
In the general formalism, we allow 7-ray photons and 
proton acceleration to occur at different locations. We 
therefore introduce a more general parameter / 7 / p (Eq(5|) 
to denote the ratio between photon luminosity and non- 
thermal proton luminosity, which reduces to f e in any 
model that invokes the same site for photon production 
and proton acceleration (e.g. the internal shock model). 
(3) We generalize the previous low optical-depth treat- 
ment to also include very high optical-depth regime by 
invoking a more general e„.i (Eq(5]). (4) We introduce 
another factor f p (EqfTU]) that represents the fraction of 
energy in those protons that can most efficiently produce 
neutrinos via the photo- pion process [6( . 

The general scheme for calculating neutrino flux is as 
follows: For an observed "Band" -function photon flux 
spectrum 
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and p is the proton spectral index defined by 
N(Ep)dE p oc E~ p dE p . The indices a v and /3 U are de- 
rived by assuming that the neutrino flux is proportional 



to the P7 optical depth t P7 , which is roughly valid when 
t P7 < 3. In this case, the first break 



7.3 x 10 5 GeV (I + z)~ 



r2.5 e 7 ,McV (^) 



is defined by the break in the photon spectrum. For the 
case of high pj optical depth, the neutrino flux no longer 
significantly increases as r p7 > 3 (essentially all the pro- 
ton energy is given to pions). The neutrino spectrum may 
be still approximately delineated as the above broken 
power law form, with e^i smaller by a factor (/tt/3)^ 7-1 . 
In general, one can write 



e„,i =e°min(l,(/j3) 1 -^), 



(6) 



where 
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is the P7 optical depth for protons with energy that in- 
teract with photons at the peak energy E? = e 7 , \' (E?) 
is the comoving proton mean free path for p^j interaction 
at EP, and AR' is the comoving width of causally con- 
nected region of the jet. The parameter f„ is also the ra- 
tio between comoving dynamical time and proton (with 
energy E^) cooling time as a result of pion production 
131). The parameter R denotes the distance of proton 
acceleration site (rather than the photon emission site if 
the two sites are different) from the central engine. The 
second break energy in the neutrino spectrum, 

e„, 2 = 3.4 x 10 8 GeV (I + z)' 1 e^ 1 ' 2 L^T^Ru (8) 

is defined by the 7r + synchrotron cooling effect, above 
which the newly produced 7r + lose energy in a time scale 
shorter than the pion decay time scale. Here e B is the 
fraction of dissipated jet energy in magnetic fields, and 
L w is the luminosity of the dissipated wind. We further 
define 
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and 



the observed neutrino number spectrum can be expressed 
as 

F (E ) — dN{E " ] 
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where E p ^ & E v $, are proton energies corresponding to 
e„,i and e u ,2, respectively (EqJSJ), and E p ^ max and E p ^ min 
are the maximum and minimum proton energy. One can 
then normalize the neutrino spectrum with the total pho- 
ton fluence, i.e. 

/ dE v E v F v (E u ) = --^[l-(l-<x P ^n>) f *} 

JO 8 J-y/p 

/•lOMeV 

x / dEyEjFj(Ej), (II) 

Jl kcV 
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where < Xp-s-tt >— 0.2 is the average fraction of en- 
ergy transferred from protons to pions, and the factor 
[1 — (1— < Xp^ir >K*] is the fraction of proton energy 
that goes to pion production, which exceeds ~ 50% when 
/tt > 3 and quickly approaches 100% at larger values. 
The coefficient 1 /8 is the product of 1 /4 (4 leptons share 
the energy of one 7r + ) and 1 /2 (on average roughly half of 
P7 interactions go to the 7r + channel when all the 7r + pro- 
cesses besides A + resonance, e.g. direct-pion production, 
and multiple pion production, are taken into account). 

777. Model-dependent neutrino flux. Below we apply 
the general formalism to different models. 

(1) Internal shock (IS) model: Both photon emission 
and proton acceleration occur at i?is ~ 10 13 — 10 14 
cm (Eql3l). One can make two simplifications. First, 
f-y/p — fe can be adopted; second, L w — £ 7 /e e can be 
adopted, where e e is the fraction of jet energy in elec- 
trons, and fast cooling is assumed. The general formal- 
ism is then reduced to the IceCube formalism |4| , except 
the additional f p correction factor and the modification 
of e„,i (EqEJ. 

(2) Dissipative photosphere (ph) model: According 
to this model, the prompt GRB spectrum is formed 
near the Thomson scattering photosphere R p h — 3.7 x 
10 11 cm £10,521^5 [16]. In order to account for the ob- 
served non-thermal photon spectrum, it is required that 
significant energy dissipation and particle acceleration 
occur at moderate optical depths [9l-lllj. The heating 
processes include small-radius internal shocks (those with 
very short variability time scales 8t <C St m in, so that 
the internal shock radius is smaller than Rph), neutron- 
proton collisional heating, or magnetic dissipation in a 
magnetized wind with a striped wind geometry. For the 
neutron-proton collisional heating scenario [lOj , small- 
radius internal shocks are needed to induce collision. In 
these shocks, significant neutrino emission is possible if 
the synchrotron-cooling-limited maximum proton energy, 

E p , max = 1.3 x 10 19 eV r^i^X^V 74 ^ 172 (where 
£ is a parameter of unity to delineate the uncertainty 
of comoving acceleration time scale with respect to the 
particle gyration radius over speed of light) exceeds the 
critical proton energy demanded by the A-resonance con- 
dition EqJTJ One can see that the condition is satisfied 
for typical parameters. For the magnetic-dissipation pho- 
tosphere model [ll|, h has been argued that the recon- 
nection process is able to accelerate protons to UHECRs 
fl9j . So for all these dissipative photosphere models, it is 
likely that both photons and cosmic rays are generated at 
Rph and significant neutrino emission is produced. Com- 
pared with the IS model, the main correction is that f n 
increases by a factor of R\a,/R p h, so that neutrino pro- 
duction is enhanced. 

One way to lower the neutrino flux for the dissipative 
photosphere model is to assume that protons are not 
accelerated to high energies to satisfy the requirement 
of Eq|T] for pion production. Particle-in-cell simulations 
show that particle acceleration is inefficient in magne- 
tized relativistic shocks [llj]. However, internal shocks 



are usually mildly-relativistic. Protons can be more ef- 
ficiently accelerated than electrons, and it is possible to 
have proton accelerated in internal shocks even if the 
magnetization parameter (the Poynting-flux-to-matter- 
flux ratio) a is as high as 0.1 [l8j]. Moreover, for a striped 
magnetic wind, internal shocks tend to facilitate recon- 
nection [20j | . which would accelerate cosmic rays to high 
energies. It seems unlikely that proton acceleration is 
completely suppressed at a dissipative photosphere. 

(3) Photosphere + internal shock (ph+IS) model: For 
any efficient dissipative photosphere model, a < 1 is ex- 
pected at the photosphere (otherwise the photosphere 
luminosity is suppressed by a factor (1 + <r)). Internal 
shocks would in any case develop at R\s > R p h, where 
a <C 1 and cosmic rays are accelerated. Even if pho- 
ton emissioon at the internal shocks may be inefficient, 
photons emitted from the photosphere would in any case 
pass through the internal shock region and interact with 
the energetic protons there to produce neutrinos. Due to 
dissipation, on average, the Lorentz factor in the internal 
shocks is expected to be somewhat smaller than that in 
the photosphere. The comoving photon number density 
in the internal shock region would be somewhat higher 
at i?is in the ph+IS model than in the IS model. This 
tends to increase the neutrino flux with respect to the IS 
model. The angular distribution of 7-ray photons in the 
proton shell comoving frame is expected to be confined 
to a narrow angle about the radial direction (because 
these photons originated at a much smaller radius), as 
opposed to isotropic distribution for the IS model. How- 
ever, since protons are roughly isotropic in the comov- 
ing frame, the resulting efficiency of neutrino production 
would not be reduced [17] . Finally, the parameter / 7 / p 
can be larger than f e , since photons are likely generated 
more efficiently in photosphere than in internal shocks. 
Considering all these factors, we expect that the pre- 
dicted neutrino flux level in the ph+IS model is roughly 
the same (within a factor of a few) as that in the IS 
model. 

(4) The internal collision-induced magnetic reconnec- 
tion and turbulence (ICMART) model and other large- 
radius magnetic dissipation models: The ICMART model 
[l2| invokes a highly magnetized outflow, which remains 
un-dissipated up to a radius Rjcmart > Ris- Emis- 
sion from the photosphere and internal shocks is greatly 
suppressed. Internal shocks help to destroy the ordered 
magnetic fields, and a strong run-away magnetic dissipa- 
tion process occurs at a larger radius -Ricmart ~ 10 15 
cm. The f„ parameter is therefore smaller by a factor 
-Ricmart/R^is ~ (10 — 100). Since the photon emission 
and proton acceleration are in the same region, one can 
take / 7 / p = f e and L w = L 1 /e e . This model therefore 
predicts a much lower neutrino flux than the IS model. 
The same applies to other magnetic dissipation models 
that invoke a large emission radius e.g. [13j . 

We calculate the neutrino flux of a typical GRB in dif- 
ferent models in Fig.l. We first calculate the neutrino 
flux within the IS model. We apply Eq|31 and adopt the 
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FIG. 1. The predicted neutrino flux for a typical GRB in 
three GRB prompt emission models: "ph": dissipative pho- 
tosphere model; "IS": internal shock model; "ICMART": 
internal-collision-induced magnetic reconnection and turbu- 
lence model. Model parameters: £7,52 = 1, 8t = 0.1 s, 
e 7 ,M C v = 0.2, a-y = 1, Ay = 2, p = 2, z = 1, T = 250, 
f e = 0.1, e B /e e = 1, i?icMART = 10 15 cm. 

following values for the measuable parameters: £7,52 = 1, 
St = 0.1 s, e 7i McV = 0-2, a 7 = 1, /3 7 = 2, p = 2, and 
z = 1. For the IS model, one can write = f e and 
L w ,52 — ^7,52/ee- We therefore have three free parame- 
ters r, f e and e B /e e . It is reasonable to take f e ~ 0.1 and 
e B /e e ~ 1. The largest uncertainty comes from T. The 
predicted neutrino flux is sensitive to T (e.g. / ff oc r~ 4 in 
the IS model). T can be constrained with several differ- 
ent methods [UHlHl, all with a lot of uncertainties. The 
IceCube Collaboration has used a "benchmark" value 
T = 300. The available data now show that T is not 
universal, and is correlated to L 1 [H, [2f| . According to 
these measurements the average value of T is somewhat 
smaller than 300, which increases the expected neutrino 



flux significantly (due to the strong T-dependence). We 
use the latest correlation [25J 

T ~ 250L°; 3 5 ° 2 , (12) 

to estimate T, which gives T ~ 250 for the example GRB. 

Also shown in Fig.l are the predicted neutrino fluxes 
for the dissipative photosphere model and the ICMART 
model (with -Ricmart ~ 10 15 cm). All other parameters 
remain the same. One can see that the ph model predicts 
a higher flux level with lower neutrino energies than the 
IS model, both favoring detections. The ICMART model 
predicts a much lower flux level than the IS model. 

IV. Current status and future prospects. The contin- 
ued search for neutrino signals from GRBs by the Ice- 
Cube Collaboration is starting to pose meaningful con- 
straints on GRB models. With the current limit, the IS 
model with the inclusion of a proper f p correction and 
r — L correlation just starts to barely violate the ob- 
servational constraint Q. The dissipative photosphere 
(ph) models are already disfavored, unless an unknown 
mechanism suppresses proton acceleration in the photo- 
sphere region. The ICMART model and other large-scale 
magnetic dissipation models are entirely consistent with 
the data. Thanks to the low neutrino background in 
the interested energy range, the upper limit would go 
down linearly with time. In a few more years, if high 
energy neutrinos are still not detected from GRBs, most 
matter-dominated GRB prompt emission models would 
be disfavored. 

It is possible that different GRBs may have different 
degrees of magnetization, so that different models may 
apply in different bursts [26| ■ A more sophisticated anal- 
ysis should invoke proper weights to different models to 
calculate the total neutrino flux. The upper limits would 
then give constraints on the weights of those models that 
invoke small emission radii. 
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